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Abstract. Wc illustrate some of the difficulties that may be encountered when 
computing photodissociation and radiative association cross sections from the same 
time-dependent approach based on wavepacket propagation. The total and partial 
photodissociation cross sections from the thirty-three vibrational levels of the h 
state of HeH"*" towards the 9 other and 6 ''ll n = 2, 3 higher-lying electronic states 
are calculated, using the autocorrelation method introduced by Heller (1978) and the 
method based on the asymptotic behaviour of wavepackets introduced by Balint-Kurti 
et al. (1990). The corresponding radiative association cross sections are extracted 
from the same calculations, and the photodissociation and radiative association rate 
constants are determined. 
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1. Introduction 

This work aims at illustrating various problems that may arise when one wants 
to compute accurate photodissociation cross sections and radiative association cross 
sections using the same time-dependent wavepacket propagation. The high sensitivity 
of the radiative association cross sections at low energy is discussed. 

While the computation of the photodissociation cross section from the lowest 
rovibrational level of the initial state is usually straightforward, the computation of 
individual cross sections for each initial rovibrational level tends not only to be time- 
consuming (as a full propagation is necessary for each v", J" v', J' transition) but also 
to suffer from additional difficulties, mostly stemming from the wider spatial distribution 
of the initial wavepacket and from the appearance of centrifugal barriers in the excited 
states potentials for high values of J'. 

We illustrate some of these difficulties through a study of the photodissociation of 
the b state of HeH+. Although HeH+ is one of the simplest closed-shell diatomic 
ions in its fundamental state, its nuclear dynamics in excited states is far from trivial 
as they are strongly coupled by non-adiabatic interactions. In the last decades, this 
cation has been the subject of numerous theoretical studies motivated by its relative 
simplicity, which allows for high-accuracy computations (Stanke et al. 2006, Miyake 
et al. 2011) ; its potential role as a diagnostic tool in plasma physics (Rosmej 
et al. 2006, Lorcau et al. 2010c) and its importance in astrochemistry (Roberge & 
Dalgarno 1982, Dalgarno 2005). It is indeed predicted to be the ffist molecular species to 
have formed in the universe, by radiative association of H+ and He (Galli & Palla 1998), 
and to be abundant in certain planetary nebulae (Cecchi-PesteUini & Dalgarno 1993) 
and Helium-rich white dwarfs (Harris et al. 2004). Every attempt of extra-terrestrial 
observation of this cation so far has however proven inconclusive at best. The short 
radiative lifetime (r ~ 10~^ s (Chibisov et al. 1996)) of its b state indicates that it 
is likely to quickly deexcitate towards the mctastable a ^E"*" state, which could play an 
important role in the astrochemistry of HeH+ because of its long radiative lifetime (r 
= 149 s for its lowest vibrational level (Loreau et al. 2010b, Loreau et al. 2013)). 

The b state exhibits properties which complicate the dynamics of its 

photodissociation. It is close in energy to a large number of other states, both bound 
and dissociative, which are strongly coupled by non-adiabatic interactions. Its large 
number of rovibrational levels also make it a good example of the importance of taking 
the vibrational dependence of the photodissociation cross sections into account. 

2. Moleculair data 

All dynamical calculations were performed using the potential energy curves, non- 
adiabatic radial couplings and dipole moments computed by Loreau et al. (2010) at 
the CASSCF and CI level using the MOLPRO quantum chemistry package (Werner 
et al. 2006). To ensure a correct description of the excited HeH+ triplet states, these 
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Figure 1. Diabatic (left) and adiabatic (right) potential energy curves of the n = 1-3 
(blue, black) and '^11 (red) states of HeH+. The b ■^S]+ state is shown in black. 

calculations were performed using the aug-cc-pV5Z basis set supplemented with one 
contracted Gaussian function per orbital per atom up to n = 4, where n is the largest 
principal quantum number of the atomic fragments. All S and 11 triplet states up to n = 
3 were included in our dynamical calculations, i.e. 11 as well as 6 ^11 states (shown 
in Fig. [T]), as photodissociation perpendicular to the laser polarization was shown to 
play an important role in the case of the fundamental state (Sodoga et al. 2009). 

Among the non-adiabatic radial couplings elements Fmm' = {Cm\dR\Cm) ^ '^^^J the 
couplings between successive states Fm,m±i were taken into account in the determination 
of the adiabatic-to-diabatic transformation matrix D{K), like in previous studies on this 
system (Sodoga et al. 2009, Loreau et al. 2010c). The equation 

dRD{R) + F{R) ■ D{R) = (1) 

was solved by continuity with the initial condition D{oo) = I, in order to ensure that the 
adiabatic and diabatic representations are identical at large intemuclear distance. The 
diabatic potential energy curves are the diagonal elements of the matrix = D~^-U^-D, 
where If^ is the matrix of the electronic Hamiltonian H'^^ in the adiabatic representation. 

As the photodissociation cross sections are highly dependent on the vibrational 
wavefunction of the initial state, we computed the wavefunctions for the b ^S"*" state 
using its potential energy curve obtained at the CASSCF and CI levels using the aug- 
CC-PV6Z basis set (instead of aug-cc-PV5Z) supplemented with the orbitals optimized 
by Loreau et al. (2010). We then resolved the vibrational nuclear equation in a basis 
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of B-Splines to obtain the energies and wavefunctions of its vibrational levels (Bachau 
et al. 2001, Loreau et al. 2010b, Cauet et al. 2012). The problem was solved for all 
possible values of J by adding the corresponding centrifugal term to the potential. 
The h ^T^' state was found to support 33 bound vibrational levels and a total of 1307 
rotational levels. 



3. Photodissociation cross sections 



The cross sections were extracted from time-dependent wavepacket propagations. In 
this approach, the initial wavepackets 0aJvv'j"(^ = 0) corresponding to the transition 
from the initial state towards the excited states m of symmetry A' are obtained by 
multiplying the nuclear wavefunction of the initial state Xa"v"J"{R) by the appropriate 
transition dipole moments ^oA",mA'{R)- 

4'0/^'v"J"{R'^t = 0) = fioA",mA'{R)XA"v"J"{R) (2) 

The wavepackets are then propagated in time, using the split-operator method 
in the present case (Feit et al. 1982, Alvarellos & Metiu 1988, Baloitcha et al. 2001). 
The total photodissociation cross section o'l°^,fj„_^/^,j,{E) can then be extracted from the 
autocorrelation function 

Col"J"(^) = ^{<PoJV'v"J"it = ^)\^OIV'v"J"it)) (3) 
m 

since it is given by (Heller 1978): 

yAE)=AnaalERe Cq^J/' (t)e'(^0''"^"+^)*/^(it, (4) 

JO 



total 

Ov"J"^A' 



where E is the photon energy, is the Bohr radius and Eq^hjh is the energy of 
the considered rovibrational level of the initial state. This method only yields the 
total photodissociation cross section, i.e. the sum of the partial, channel-specific cross 
sections: 

(7o°''J"^-A'J'(-^) = X! '^Oi'"J"-s>mA'J'(-^) (^) 
m 

The total cross section characterizes the depopulation of the initial state caused 
by the absorption of a photon rather than the dissociation of the molecule into specific 
fragments. It is therefore unsuitable for the computation of branching ratios or of 
radiative association cross sections from a specific excited channel towards the initial 
state (see Section |4]). 

The cross sections exhibit predissociation resonances caused by the non-adiabatic 
crossing between a bound state and a lower dissociative channel: trapped wavepackets 
move back and forth in the potential energy well of the bound state until they non- 
adiabatically cross to the dissociative channel, leading to the appearance of resonances 
in its spectrum (Balakrishnan et al. 1999). 
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We chose to neglect them completely in the present calculations for several reasons. 
Their exact determination would have required impractically long propagation times, 
as a non-negligible fraction of the initial wavepacket was still found in the potential 
energy wells after propagation times as high as 50 x 10^ a.u. of time (1.2 ns). 
Moreover, taking these resonances into account may lack physical sense if one wants 
to extract photodissociation and radiative association cross sections from the same 
wavepackets propagations, as discussed in Section |4j Furthermore, including very 
long-lived resonances in our photodissociation cross sections may be unrealistic as 
electronic transitions caused by spontaneous radiative deexcitation, collisional processes 
or photoexcitation are likely to take place on shorter timescales. Finally, the impact of 
Fano-type resonances on the photodissociation rate constants is probably small (see 
Section |5]). Such resonances may be more adequately described using partitioning 
techniques (Desouter-Lecomte & Lievin 1997) or the time-independent formalism (van 
Dishoeck & Visser 2011, Desouter-Lecomte et al. 1995). 

Predissociation resonances can easily be suppressed when computing the total 
photodissociation cross section through the autocorrelation method for low values of 
v. their appearance occurs on a longer timescale than direct photodissociation, as 
illustrated by the autocorrelation function in the case of the photodissociation of the 
V = 5 level (Fig. [2| black curve). The departure of the wavepackets from the Franck- 
Condon region and the return of their trapped fractions occur on two distinct timescales 
and, if the propagation is stopped before their return, the computed cross section shows 
no sign of resonances (Heller 1978) and exhibits instead a smooth continuation of the 
cross section below the threshold energy. This can be understood intuitively as the 
autocorrelation function then contains no information about the trapped fraction of 
the wavepackets and behaves as if the totality of the wavepackets could directly reach 
the asymptotic region (Balakrishnan et al. 1999). This smooth envelope gradually 
disappears for longer propagation times as the Feshbach resonances appear in its 
stead. While a resonance-free total cross section can easily be obtained for the first 
few vibrational levels by using short propagation times, it ceases to be applicable 
for higher vibrational levels as the Franck-Condon region becomes wider, keeping the 
autocorrelation function from ever falling to zero (Fig. [2| red curve). 

Partial photodissociation cross sections can be computed using the method 
introduced by Balint-Kurti et al. (1990). In that approach, the cross sections are given 
by the Fourier transform of the component of the wavepackets reaching an internuclear 
distance -Roo located in the asymptotic region: 

partial { JP\ — ^L^^l'^^^^^ JP\ A^^" J" { JP'\\'^ 



/i 



with 



A'S/;iE) = ^r 0^^;/;(i?oo, t)e^(^o."."+i^)*Arft (7) 

V27r "'0 



where fj, is the reduced mass of the molecule, a is the fine structure constant and 
km = J'^fJ'{EQy//j" + E — -E^ymptotic) is the magnitude of the wave number in the 
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Figure 2. Autocorrelation functions for the pliotodissociation of tire v = 5 (black) 
and w = 15 (red) levels of the b ^S"*" state towards higher excited "^S"'" states. In 
the case of w — 5, two different time scales can clearly be observed (with the direct 
photodissociation contribution appearing before 2000 a.u. of time), unlike in the case 
of t; = 15. 



considered channel m. 

The propagations were performed on a spatial grid of 2^^ points that spans from an 
internuclear distance i?mm = 0.1 a.u. up to -Rmax = 200 a.u. The use of a long spatial 
grid was made necessary by the components of the wavefunctions at large internuclear 
distance for the highest initial vibrational levels. The Fourier transform was performed 
at Roo = 175 a.u. and we placed a quadratic optical potential starting ai R = 180 
a.u. in order to avoid reflections of the wavepackets at the edge of the grid. Since this 
method requires the wavepackets to reach the asymptotic region, it necessitates longer 
propagation times than the computation of the total photodissociation cross section 
through their autocorrelation function. In the present work, propagations up to 5 x 10''^ 
atomic units of time (1.2 ns) were performed with a time step of 1 a.u. for all 33 initial 
vibrational levels. Tests with shorter time steps were carried out to ensure accuracy. 

Unhke the autocorrelation-based approach, this method yields state-specific 
information. Moreover, it only takes into account the parts of the wavepackets that reach 
the asymptotic region of the dissociation channels and ignores their trapped fractions 
as long as they stay in the potential energy wells. As a consequence of this, the partial 



Time-dependent determination of photodissociation and radiative association cross sections! 




1.70 1.75 1.80 1.85 1.90 

Photon energy E^^(eV) 

Figure 3. The Gibbs phenomenon leads to the appearance of discontinuities near 
the threshold energy in the partial cross section from the v = 2 level of the h ^Ti^ 
state towards the H+ + He(ls2p '^P°) fragments (full line, black), while the total cross 
section as obtained through the Fourier transform of the autocorrelation function for 
short propagation times is perfectly smooth (dashed line). Thanks to this, the correct 
form of the partial cross section can be deduced (full line, red). 

cross sections obtained with this method often start abruptly at the threshold energy and 
are thus affected by the Gibbs phenomenon (Jerri 1998), which causes ringing artifacts 
around jump discontinuities in Fourier transforms, as illustrated in Fig. |3| 

Such Gibbs oscillations only appear if the cross section is non-negligible at the 
threshold energy. The reflection principle allows for an intuitive understanding of the 
conditions under which this will happen: the dissociation channel has to be a bound 
electronic state, and its bound part needs to be in a region where the wavefunction of 
the initial rovibrational state is non-negligible (Fig. |4| cases (c) and (d)). This is more 
likely to be the case for excited vibrational levels because of the wider spread of the 
vibrational wavefunction (Fig. |4| compare cases (b) and (d)). Despite this, the cross 
section may not start abruptly at the threshold energy for certain specific values of v" 
as one of the nodes of the wavefunction may coincide, through the reflection principle, 
with the asymptotic energy of the dissociation channel (Fig. |4| case (e)). Whether 
or not Gibbs oscillations appear thus depends on several factors: the topology of the 
potential energy curves, their relative positions and the initial vibrational level. 
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Figure 4. Schematic illustration of various cases where the partial photodissociation 
cross section a does not start abruptly at the threshold energy (a, b, e) and where 
it does (c, d), depending on the shapes and relative position of the potential energy 
curves and on the initial vibrational level. 



Although various filtering methods may be used to reduce these ringing artifacts 
(Jerri 1998), it is generally impossible to suppress them completely. It is however 
possible to do so in the present case using the total photodissociation cross section 
calculated from the autocorrelation function for short propagation times (Fig. |5|: since 
it (incorrectly) describes predissociation as a smooth envelope, no jump discontinuity 
occurs at the photodissociation treshold energy. It is therefore unaffected by the Gibbs 
phenomenon (Fig. |3]). As the sum of the partial cross sections must be equal to the total 
cross section in the energy ranges where it contains no contribution from the incorrectly- 
described predissociation, the total cross section obtained for short propagation times 
can thus be used as a reference to suppress the ringing artifacts in a given partial cross 
section for a channel m'^: 

^Ov"J"^mfk'J'{^)=^o'v'^]"^k'J'{^)~ X! '^Oi>"J"-s>m'A'J'(-^) {^) 

However, this is only possible if the ringing artifacts do not occur in the same 
range of energy for two different exit channels, as their respective contributions to the 
total cross sections cannot be separated in that case. By combining the results of our 
computations of the partial and total photodissociation cross sections, we were able to 
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Figure 5. Illustration of the reflection principle in the photodissociation cross section 
of the u = 5 level of the h state (black) towards the H+ + He(ls2p '^P°) 

channel (red). The Fourier transform of the autocorrelation not only yields the 
direct photodissociation cross section (green, full line), but also a contribution from 
predissociation which incorrectly appears as a smooth envelope for short propagation 
times (green, dotted line). 



reduce or suppress the ringing artifacts near the threshold energy and to obtain the 
partial photodissociation cross sections from the 33 vibrational states of the h state 
through the 15 other n = 2-3 ^S"*" and ^11 dissociation channels. 

The impact of the vibrational excitation of the initial state on the photodissociation 
cross sections is twofold: the nodal structure of the wavefunction directly affects the 
photodissociation cross sections (compare Fig. |4] (b) and (d)), moreover, the spread of 
the wavefunction towards high values of R tends to favor transitions towards different 
excited channels. 

Partial photodissociation cross sections free from Feshbach resonances can often 
be obtained since the fraction of the wavepackets that dissociates directly usually 
reaches the asymptotic region on a much shorter timescale than the fraction that stays 
momentarily trapped in potential energy wells. This is however not always true when 
considering the photodissociation of excited vibrational levels. 

Shape resonances caused by tunneling through the centrifugal barriers are also 
expected to appear for high values of J', but the very large number of rotational levels 
of the b ^S"*" state precluded us from computing individual cross sections for each of them 
as a full wavepacket propagation would be necessary for each of the J" — )■ J' = J" ± 1 
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transitions towards other states and J" — )■ J' = J", J" ± 1 transitions towards 
states. Wc therefore chose to neglect the rotational dependence of the cross sections, as 
its impact on the cross sections was shown to be smaller than that of the vibrational 
dependence in the case of the X ^J]^ state (Sodoga et al. 2009, Miyake et al. 2011). 
However, although the envelopes of the cross sections may not vary significantly with J" 
and J', shape resonances have been shown to significantly augment photodissociation 
and radiative association rate constants in some systems, especially at low temperature 
(Antipov et al. 2009). It is common to compute the envelope of the cross section 
through time-dependent methods without regards to shape resonances as their exact 
determination would require an impractically fine energy grid and long propagation 
times. They can indeed be computed separately, e.g. through the Breit-Wigner 
formalism (Breit & Wigner 1936, Bain & Bardsley 1972). 



4. Radiative association Cross Sections 



Radiative association cross sections are most commonly obtained through time- 
independent calculations, which require the discretization of the energy continuum and 
separate calculations for each channel. By contrast, Martinazzo and Tantardini proposed 
a wavepackct-bascd method of computation of the radiative association cross sections, 
motivated by the superior scaling properties of time-dependent methods with regards to 
the size of the system, the number of channels considered and the energy range studied: 
a single time-dependent calculation yields the radiative association cross section towards 
all channel and for several energies at once (Martinazzo & Tantardini 2005). 

We chose a different time-dependent approach, which retains these advantages: 
since radiative association is the inverse process of photodissociation, the transition 
dipole matrix elements Mq^,, j„ jt — w J"\^J|'\{^ mv' J')\'^ involved in both processes 
arc the same, with the photodissociation cross section from the initial level v" , J" of the 
state being given by (Barinovs et al. 1999) 

8 TT^ ■^"^^ 

(ylv°^J^^^^rf^^^'^\Ehv) = 3 [A-Keo) ^ ^^J",J'^lv"J",mv>J'{E) (9) 

and the corresponding radiative association cross section from the channel m towards 
all v" , J" levels of the state being given by (Barinovs et al. 1999) 

64 TT^ "^"^^ n 

= y (4vreo) 5' ^^^"'^'^^'^''^''-'^'^^^^^^^ 

where Sjn^j' are the Holn-London factors and Qm is the degeneracy factor of the 
considered channel. 

The computation of both partial photodissociation and radiative association cross 
sections thus comes down to the determination of the same transition dipole matrix 
element Mq^h jh .^^i j, = \{'^ov"j"\lJ'\{'^mv'j')\'^ between a bound and a free state. 
The radiative association cross section can thus directly be obtained from the same 
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wavepacket propagation as the corresponding photodissociation cross section or from 
the photodissociation cross section itself. By comparing Eq. ^ and (10), it is indeed 



seen that the radiative association cross section a^^l'^n^},; from an excited channel m 
towards the v", J" level of a lower state is related to the cross section characterizing the 
photodissociation from that v", J" level of the lower stat^ towards the channel m by 
the relation (Puy et al. 2007): 



2 



Rad.Stab. fei^ ^Photodiss . .partial 

'^mv'J'^Ov"J" — ,,„2Z? '^Ov"J"^mv'J' \^^) 

While the photodissociation cross sections are expressed as a function of the incident 
photon energy E^^, the radiative association cross sections are a function of the relative 
kinetic energy E/^ of the two colliding fragments. The two energy scales differ simply by 
the photodissociation threshold energy £'thresh,?;"j"— i-6. the energy difference between 
the initial level v"J" and the asymptotic energy of a given fragmentation channel m: 

Ek,v"J"^m = Ehu — EthTesh,v"J"^m (12) 



The presence of the collision energy E^ in the denominator of Eq. [TT] implies that 
radiative association cross sections at low energy, and therefore the rate constants at low 
temperature, are highly sensitive to the value of the corresponding photodissociation 
cross sections near the threshold energy. The suppression of ringing artifacts is thus 
important to the determination of accurate radiative association cross sections and in 
cases where the corresponding photodissociation cross sections are non-negligible near 
the treshold energy, as illustrated in Fig. |4| the radiative association cross sections will 
be larger for low collision energies. 

Although computing a radiative association cross section through the propagation 
of dissociative wavepackets may seem counterintuitive, it presents some advantages: it 
circumvents the problem of the choice of the shape of the initial wavepacket in collisional 
problems and it allows the determination of the cross section on its whole range of 
energy through a single propagation, whereas several propagations are necessary to 
cover a wide range of energy in a time-dependent collisional approach (Martinazzo & 
Tantardini 2005). Note, however, that the dissociation approach used here describes 
the radiative association from several collision channels towards a single bound state, 
whereas the usual collisional approach describes the radiative association from a single 
collision channel towards several lower states at once. Our approach is more time- 
efficient if one wants to study the formation of a molecule in a specific state by radiative 
association from several collision channels at once. The usual collisional approach is 
however preferable if one wants to study the different bound electronic states that can 
be formed by radiative association from a specific initial channel. 

Radiative association can occur towards any of the rovibrational level of the inferior 
state, its cross section a^^^^^^' is thus obtained by summing the contributions of 

I For clarity, we use the following convention: u", J" denote the rovibrational level of the lower 
electronic state (i.e. the initial state in photodissociation but the final state in radiative association) 
while v', J' denote the rovibrational level of the upper state. 
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Figure 6. Cross section for the radiative association of He+ (Is) and H(ls) towards 
the fundamental X state via the A ^S"*" channel, with (full line) and without 
(dotted lines) taking its vibrational dependence into account. 

all the rovibrational levels <ym1y'P-^v"j"- Since the rotational dependence of the cross 
sections was neglected in the present study, we approximated the radiative association 
cross sections by summing the cross sections for all vibrational levels multiplied by the 
corresponding number of rotational levels. 

We tested our method on the well-studied case of the photodissociation and 
radiative association between the X ^S"^ and A states of HeH"*" (Zygelman & 
Dalgarno 1990, Kraemer et al. 1995, Miyake et al. 2011). In this specific case, we 
found the neglect of the vibrational dependence in photodissociation to be a very bad 
approximation if one uses them to compute the corresponding radiative association 
cross sections (Fig. |6]): since the photodissociation cross section for v" = is very small 
near the threshold energy for the A ^S"*" state (~ 13 eV, while the maximum of the 
cross section is located around 25 eV) (Roberge & Dalgarno 1982), its contribution to 
the radiative association cross section at low energy is particularly small. The largest 
contributions to the radiative association cross section actually come from the transitions 
involving the v" = 7 and v" = 8 levels of the X state, as the corresponding 
photodissociation cross sections are large near the threshold energy. 

While not as dramatic, the vibrational dependence also has a significant impact in 
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Figure 7. Cross sections for tlie radiative association of He(ls3s '^S) + H+ (black 
curves), He(ls3p ^P°) + H+ (red curves), and Re{ls3d '^S) + H+ (blue curves) 
towards the b ^Tj~^ state of HeH+ with (full lines) and without (dotted lines) taking 
the vibrational dependence of the cross sections into account. 

the case of the radiative association of the b state (Fig. [?]), modifying the cross 
sections by several orders of magnitude. 

Our radiative association cross sections were obtained by suppressing the Feshbach 
resonances, as the non-adiabatic transitions they are caused by are unhkely to occur 
in the course of a colhsion. This can be understood intuitively in the time-dependent 
approach we use: during photodissociation, wavepackets trapped in a bound excited 
state keep oscillating in its potential energy well until they non-adiabatically cross to 
the continuum of a lower-lying state (provided no other phenomenon takes place, such 
as collisions or radiative decay). The non-radiative transition thus has repeated chances 
to occur, whereas it only has a small probability of taking place during a coUisional 
process such as radiative association. 

5. Rate constants 

The rate constants for the photodissociation of the initial state caused by photons 
emitted by a blackbody of temperature and of radius R^, situated at a distance R, 
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is given by (Roberge & Dalgarno 1982): 

^P.tot / n^\TP2 



P 47r (R.\' j al^%{E,^,T)El 



where aQ^°l,{Ehi,,T) is the total photodissociation cross section characterizing the 
destruction of the state assuming a Maxwell-Boltzmann distribution of population 
in the rovibrational levels: 

^o-.A' [Eh.. T) - ^^,,,„(2J" + 1) e-^^",^"A-^ ^^^^ 

Note that T stands for the temperature of the Maxwell-Boltzmann distribution 
of population in the initial state while T* denotes the temperature of the blackbody 
emitting the photons responsible for the photodissociation process. 

As in Roberge and Dalgarno's study of the photodissociation of the X^T,'^ state of 
HeH+, we chose a ratio {E^/Kf of 10~^^, which corresponds to the physical conditions 
met in planetary nebulae such as NGC7027. As previously mentioned, the cross sections 
were assumed to be independent of J" (i.e. (Jv,,}" = o'v,j"=o for all J"). The impact of 
the Fano-type resonances on the photodissociation rate constants was not assessed, but 
it is expected to be small, as their asymmetric profile would likely limit their net result 
upon integration of the cross sections. 

As in all calculations of this kind, it should be noted that the photodissociation 
rate constants are likely to be underestimated at high blackbody temperature due to 
the absence of excited states beyond n = 3 in our computational basis. 

Likewise, the radiative association rate constants were computed assuming a 
Maxwell-Boltzmann distribution of the relative kinetic energy of the colliding fragments: 



k 



RS 



I ^^f^^^''j"E,e-^'/''^^dE (15) 

Both photodissociation and radiative association rate constants vary strongly with 
the initial/final vibrational level considered. Qualitatively similar results are to be 
expected for all systems in which the photodissociation cross sections of the first 
vibrational levels differ significantly, which depends both on the shapes and on the 
relative positions of the potential energy curves. Whether neglecting the vibrational 
dependence of the cross sections is a valid approximation or not is therefore expected 
to vary from one system to another. 



6. Conclusions 



We have illustrated the main difficulties that can arise when calculating photodisso- 
ciation cross sections using time dependent methods, particularly when one wants to 
deduce radiative association cross sections from the same calculations. We used the 
photodissociation of the b of HeH"*" as an example as this state exhibits several 
peculiarities, such as a large number of vibrational levels and a relative proximity to 
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Figure 8. Rate constants for photodissociation of the h state as a function of 
the temperature of the blackbody, for different Maxwell-Boltzmann distribution 
temperatures. 

several excited electronic states strongly coupled through non-adiabatic interactions. 
The presence of bound states among those was shown to further comphcate matters. 
Firstly, it leads to the occurrence of Feshbach resonances, which we chose to neglect 
here since their exact determination through time-dependent approaches is impractical 
and since the corresponding non-adiabatic transitions are unlikely to occur during ra- 
diative association. Secondly, it causes Gibbs oscillations to appear near the threshold 
energy in partial photodissociation cross sections, which significantly affect the radiative 
association cross sections at low collision energy. 

By combining two methods of computation of the photodissociation cross section, 
the partial photodissociation cross sections from the thirty-three vibrational levels of 
the h ^E"*" state of HeH"^ towards the 9 other ^Ti^ and the 6 ^11 = 2 and 3 higher- 
lying electronic states have been calculated with minimal contributions from the Gibbs 
phenomenon. Thanks to judicious choices of the propagation time, photodissociation 
cross sections free from Feshbach resonances were obtained in order to compute the 
corresponding radiative association cross sections. We have shown that this approach 
may be more time-efficient than the coUisional approach to describe the radiative 
association towards a specific state for several collisional channels at once. The 
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photodissociation and radiative association rate constants were computed and were 
shown to be strongly affected by the vibrational dependence of the cross sections. 
Although it is expected to be smaller, the role of the rotational dependence of the 
cross sections could be significant due to the contributions of shape resonances. 
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